Inter-dimensional Degeneracies in van der Waals Clusters and 
Quantum Monte Carlo Computation of Rovibrational States 



M. P. Nightingale 

. Department of Physics, University of Rhode Island, Kingston RI 02881, USA 

o 
o . 

. Mervlyn Moodley 

C . 

School of Physics, Howard College, University of KwaZulu- Natal, Durban, ^0^1, South Africa 



(N 

> 



a: 



(Dated: February 1, 2008) 

Quantum Monte Carlo estimates of the spectrum of rotationally invariant states 
of noble gas clusters suggest inter-dimensional degeneracy in N — 1 and N + l spacial 



dimensions. We derive this property by mapping the Schrodinger eigenvalue problem 
onto an eigenvalue equation in which D appears as a continuous variable. We discuss 
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implications for quantum Monte Carlo and dimensional scaling methods. 
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I. INTRODUCTION 

One of the advantages of Monte Carlo methods is that they scale well with the number 



of degrees of freedom of a physical system. In this paper we consider van der Waals clusters 
consisting of N bosonic Lennard- Jones atoms in D spatial dimensions. We treat the atoms 
as "elementary" particles without internal degrees of freedom, so that in total we deal with 
clusters with ND quantum mechanical degrees of freedom. We are mainly interested in the 
energy spectra of these clusters. 

Quantum Monte Carlo computations can be made much more efficient by the use of 
optimized trial wave functions, as is well known. With currently available methods, as a 
matter of fact, the problem of computing rovibrational spectra with Monte Carlo methods 
is virtually intractable without good trial functions. One of the questions of interest is the 
relative importance of the quality of these trial wave functions for n-body correlations with 
n in excess of the commonly used correlations with n = 2 and n = 3. In this context, 
the idea of varying the spatial dimensionality of the system quite naturally suggests itself, 
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because particles can be more compact in higher dimensions, which suggests that correla- 
tions involving a higher number of particles might become more important as the spatial 
dimensionality increases. While we have not found clear numerical evidence to support this 
idea, 1 our computations did produce an interesting by-product, which forms the topic of 
this paper. 

Our computations showed that the energy spectra of N particles in N — 1 and N + 1 
spatial dimensions are numerically indistinguishable for states invariant under rotation and 
translation. 1 Indeed, in this paper we show that for these S'-states and for D > N — 1, the 
N particle time-independent Schrodinger equation can be transformed into an eigenvalue 
equation involving a differential operator with |(7V — 1)N independent variables and an 
effective potential in which the spatial dimension D appears as a continuously varying pa- 
rameter. This effective potential turns out to depend quadratically on D and is symmetric 
about D = N. This implies the aforementioned inter-dimensional degeneracy, in addition 
to a relationship between the wave functions for D = N — 1 and D = N + 1 dimensions. We 
note here that, as is also manifest in the frustration of the classical system for D < N — 1 
—frustration for example in the sense that not all inter-atomic distances can be equal— 
the spectrum for values of D < N — 1 cannot be obtained by analytic continuation of the 
spectrum for D > N — 1; we shall return to this in the discussion at the end of this paper. 

Inter-dimensional degeneracy was also derived recently by Gu et al, 2 by a group- 
theoretical method. To the best of our knowledge, inter-dimensional degeneracies of S'-states 
of Lennard- Jones clusters have not been observed before, with the exception of the two-body 
cluster in one and three dimensions which follows from the standard separation of variable 
solution of the Schrodinger equation. 

Even though an extension of our study to dimensions higher than the physical three 
dimensions is primarily of academic interest, the effect of spatial dimension on quantum 
systems has been studied since the early days of quantum physics. In fact, Fock 3 as early as 
1935 showed that there exists a relationship between the the hydrogen-like wave functions 
and four- dimensional hyper-spherical harmonics. 4 The hyper-spherical coordinate method 
was used in the late seventies to discover inter-dimensional degeneracies in electron systems. 
For the one-electron system a transformation was found that reveals inter-dimensional de- 
generacy between a system in D dimension and angular momentum / with the same sys- 
tem in D ± 2 dimensions and angular momentum I =p l. 5 Many-electron systems have also 
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been shown to exhibit inter-dimensional degeneracies. 6 Further references to other inter- 
dimensional studies can be found in Ref. 2. 

The structure of this paper is as follows. In Sec. II A we briefly summarize the quantum 
Monte Carlo technique for excited states developed in Refs. 7 and 8 to obtain optimized trial 
wave functions for van der Waals clusters. In this method, trial wave functions are developed 
that can be improved systematically. These trial wave functions are linear combinations of 
elementary basis functions with non-linear variational parameters. The elementary basis 
functions consist of a prefactor and exponential polynomial that is formulated in terms of 
all possible iV-body correlations. In Sec. II B we present Monte Carlo energy estimates 
obtained for selected few-body van der Waals clusters in a limited number of dimensions 
ranging from D — 1 to D — 6. In this study we consider van der Waal clusters composed 
of atoms of Kr, Ar, Ne and the hypothetical §-Ne, which has half the (dimensionless) mass 
of Ne. Kr can be considered as a semi-classical case while that of §-Ne is more quantum 
mechanical in nature. Sec. Ill is devoted to the exact derivation of dimensional degeneracy, 
with some of the results postponed to the Appendix. In the final Section IV we discuss the 
relevance of our results, in particular for dimensional scaling methods. 



II. MONTE CARLO APPROACH 



A. Optimization of ground and excited state wave functions 

We consider clusters in D dimensions consisting of N atoms with positions specified by 
the D x N matrix of Cartesian coordinates R = (rir2 . . . with 



(1) 



We shall use the following definitions 



1 u 



Tj ~ r { 



(2a) 
(2b) 



for difference vectors and their lengths. 
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For a system of iV bosonic van der Waals atoms with atomic mass \x and interacting via 
a pair potential, the dimensionless Hamiltonian is 

N 

■In) 



1 N 

^ = -^£ v * + £^)> ( 3 ) 



<=1 (hj) 

with 



V ? = Et£- (4) 



i 9x L 

Q = l " l 



and where 1/ is the dimensionless Lennard- Jones potential 



V(r) = ± ~ ^ (5) 



The inverse dimensionless mass is given by vrT x = h 2 /2s pa 2 e, which is proportional to the 
square of the de Boer parameter, 9 where — e is the minimum of the Lennard- Jones potential 
and 2^ a the corresponding inter-particle distance. 

A preliminary step in our optimization procedure is to generate a sample of configura- 
tions R CT , with a — 1, . . . , s, which are sampled from a relative probability density function 
^(Ro-) 2 . The guiding function ip g used for the computations reported in this paper is defined 
in terms of a trial function tp, which approximates the ground state. In simple cases, we 
used i\) 2 g = ?p 2 / p , where the parameter p is chosen in the range 2 < p < 3. Where necessary, 
we used a more sophisticated guiding function 10 so as to generate a sample with substantial 
overlap with all the excited states under consideration. 

The trial wave functions are linear combinations of elementary basis functions fa, each of 
which implicitly depends on non-linear variational parameters, and we use different proce- 
dures to optimize the linear and non-linear parameters. For reasons explained in detail below, 
we define the re- weighted functions ^(R) = i) g {R)~ l f3i{R) and ^(R) = i) g {R)~ 1 H For 
a complete set of elementary basis functions fa the Schrodinger equation can be written in 
the form 

n 

A'(R CT ) = ^(R^- (6) 

In practical applications, the set of functions fa is, of course, far from complete, but the nxn 
matrix £ may still be determined by solving Eq. (6) for S in a least-squares sense given the 
re-weighting just introduced. Note that Eq. (6) is exactly satisfied if the functions fa span 
an invariant subspace of the Hamiltonian H, even if they do not form a complete set; this 
provides an important zero-variance principle for the corresponding part energy spectrum. 
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It is convenient to rewrite Eq. (6) in matrix form 

B' = B£, (7) 

where B ai = fa(R a ) and B' ai = Multiplying Eq. (7) from the left by the transpose 

of B, one obtains by inversion 

£ = (B T B) _1 (B T B / ) = N _1 H, (8) 

with Nij = h{Ra)(3j{Ro) an d Hij = J2<r Pi(R<r)Pj(R<r)- As can be easily verified, Eq. (8) 
is indeed the least-squares solution of Eq. (6). Note that for an infinite sample the hermiticity 
of the Hamiltonian guarantees that H is a symmetric matrix, but this is not the case for a 
finite Monte Carlo sample. If H is symmetrized in Eq. (8), the resulting £ no longer satisfies 
the least-squares property nor the aforementioned zero- variance principle. 

The optimal linear combinations of the basis functions fa are computed by constructing 
the spectral decomposition of £: 

n 

£ l3 = Y^d\E k d) (9) 
k=i 

where dj and d\ are the components of the left and right eigenvectors of £ with eigenvalues 
E k . This yields the trial functions 

n 

^* = (io) 

i=l 

Before we continue this review of our optimization procedure, some comments should 
be made. First of all, the matrix N frequently is ill-conditioned. This problem can be 
dealt with by using a singular value decomposition of the matrix B to obtain a numerically 
regularized inverse of N. 7 ' 8 Secondly, we mention that Eq. (10) can also be derived from 
the condition that the quantum mechanical expectation value of the Hamiltonian in the 
states ip k is stationary with respect to variation of the coefficients d\. This condition yields 
a generalized eigenvalue equation involving matrices N and H, the Monte Carlo estimators 
of which are the matrices N and H introduced previously. The re-weighting defined before 
Eq. (6) was introduced so that these estimators are unbiased. 

As mentioned, the basis function fa depend implicitly on non-linear variational parame- 
ters. These are optimized iteratively and it should be kept in mind that for each choice of 
the non-linear parameters, new optimized linear parameters have to be computed. The full 
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optimization of all parameters therefore entails a linear optimization nested in a non-linear 
one. The linear optimization is a standard linear algebra problem; the optimization of the 
non-linear parameters is performed by minimizing the variance of the local energy of the 
wave function: 

= e: =1 [^(r.)-^(r ct )] 2 (n) 

where i[) k = ipg 1 ^ and ijj kl = il)~ 1 Hi[) k . 

The trial wave functions produced by this method yield estimates of the energy levels 
that are upper bounds to the exact energies, if statistical errors are negligible. To reduce 
the resulting systematic errors, the so-called variational errors, we employ these optimized 
wave functions as basis functions in a correlation function Monte Carlo calculation. 11-13 This 
reduces the variational bias in the eigenvalue estimates, but it usually increases the statistical 
errors in the estimates. In a formal sense, the reduction of variational errors obtained in 
correlation function Monte Carlo is accomplished by introducing a new and improved basis 
by means of the substitution 

Pi(R) - exp(-tH)/3i(R) = fo(R,t). (12) 

For increasing projection time t the spectral weight of more highly excited states in the 
new basis is reduced, and with it the variational error. In the limit t — > oo all states of 
the new basis collapse onto the ground state, which implies that as t increases, the overlap 
matrix of the t-dependent basis states becomes more nearly singular, which increases the 
statistical errors. In principle, the errors increase exponentially; in practice, the method 
as we currently use it, breaks down once the Monte Carlo estimate of the overlap matrix 
develops negative eigenvalues. 

We use elementary basis functions of the following general form 14 

Pi(R) = Si (R) exp(J2 ajSj(R) + A ( r -r)), (13) 

j cr<T 

where the term involving A imposes short- and long-range asymptotics; the Sj and Sj are 
bosonically symmetrized monomials. The exact structure of these basis functions is of no 
concern in this paper. A detailed description of the above mentioned method and the 
structure of the basis functions can be found in Refs. 7,8. 
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B. Numerical results in various dimensions 

In this section we present numerical results that show that the energy spectrum as a 
function of dimensionality for D > N — 1 is symmetric about D = N. We discuss results 
for Kr, Ar, Ne and the hypothetical |-Ne, which are defined respectively by the following 
inverse masses: 1.9128 x 1(T 4 (Kr) 6.9635 x 1(T 4 (Ar) 7.0920 x 10~ 3 (Ne) 1.4184 x 10~ 2 
(J-Ne). 

1. The three-body case 

Table I shows the ground state energies E\ for Kr 3 , Ar3 and 7|-Ne 3 in dimensions ranging 
from D — 1 to D — 6. We fitted the computed values for D > 2 to a parabola with its 





Kr 3 


Ar 3 


i-Ne 3 


D 


Ei 


AEi 


Ex 




Ei 


AEi 


1 


-1.872 548 547 6 


-9x10" 


-1 


-1.734 808 71 


-8x10" 


-1 


-0.895 584 


-4x10- 


-1 


2 


-2.760 461 351 5 


2xl0~ 


10 


-2.552 953 22 


-1x10" 


-9 


-1.302 484 


-7x10- 


-7 


3 


-2.760 555 278 7 


6xl0~ 


10 


-2.553 289 43 


1x10" 


-8 


-1.308 442 


9x10- 


-6 


4 


-2.760 461 351 3 


-5x10" 


11 


-2.552 953 22 


-1x10' 


-9 


-1.302 483 


-2x10" 


-6 


5 


-2.760 179 569 8 


-1x10" 


-9 


-2.551 944 61 


-2x10" 


-8 


-1.284 627 


-1x10" 


-5 


6 


-2.759 709 937 6 


5xl0~ 


10 


-2.550 263 64 


7x10- 


-9 


-1.254 901 


5x10" 


-6 



TABLE I: Ground state energies E\ (with errors in the least significant digit) and deviations from 
quadratic fits AEi for Kr 3 , Ar 3 and |-Ne 3 in dimensions D = 1 through D = 6. 

minimum at D — 2. The difference between the computed and fitted results AEi is also 
shown in the Table I. As is the case with the classical minimum of the energy, which equals 
—2.03 for D — 1 and —3 for D > 2, the quantum mechanical ground state energy at D — 1 
is nowhere near the curve. 

2. The four-body case 

Table II shows the ground state energies E\ for a four-body cluster, Ar 4 , in various 
dimensions. From these results it is evident that an inter-dimensional degeneracy exists 
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in D = 3 and D = 5 dimensions. Again, AEi represents the difference of the computed 
energies and the results obtained from a parabolic fit with its minimum at D = N = 4, this 
time excluding D = 1 and D = 2. The quantum mechanical estimates can be compared 
with the classical minima, viz., —3.07 for D — 1, —5.07 for D = 2, and —6 for D > 3. 





Ei 




1 


-2.625 622 56 


-2x10" 


-0 


2 


-4.329 517 95 


-8xl0~ 


-i 


3 


-5.118 146 05 


-2x10" 


-9 


4 


-5.118 653 84 


3x10" 


-9 


5 


-5.118 146 05 


-2x10" 


-9 


6 


-5.116 622 70 


1x10" 


-9 



TABLE II: Ground state energies (with errors in the least significant digit) and deviations from 
quadratic fits AEi for Ar4 in dimensions D = 1 through D = 6. 

3. Excited states 

Thus far we have only numerically verified that inter-dimensional degeneracies exist for 
ground state energies. Table III and IV strongly suggest that the same holds for excited 
states. The first table shows the four lowest excited state energies obtained for Ar 3 in 
D = 2, 3 and 4 dimensions. Once again, as observed for the ground state energies, these 
degeneracies exist in this three-body cluster for the D = 2 and D = 4 case. The D = 3 case 
is included in this table to indicate that the energies obtained here are different and lower 
than the other two cases. The second table shows the four lowest excited state energies 
obtained for Ar 4 in D = 3 and 5 dimensions; we denote energy levels by E ± < E 2 < • • • . 

The results in TableV illustrate the loss of accuracy that occurs for five particle clusters. 
The differences between the estimates of the energies of corresponding levels for four and six 
dimensions are due to the failure to converge of the correlation function Monte Carlo. This 
reflects the fact that our trial wave functions can in principle be systematically improved 
only for cluster sizes N < 4, because they contain fully adjustable n-body correlations with 
n < 4 only. 
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k 


D = 2 


D = 3 


D = 4 


2 


-2.249 860 2 


-2.250 185 5 


-2.249 860 


3 


-2.126 038 8 


-2.126 361 


-2.126 039 


4 


-1.996 153 


-1.996 43 


-1.996 153 


5 


-1.946 3 


-1.946 7 


-1.946 3 



TABLE III: Comparison of the excited state energies Eu (with errors in the least significant digit) 
of Ar3 in D = 2, 3 and 4 dimensions. 



k 


D = 3 


D = 5 


2 


-4.800 897 73 


-4.800 897 75 


3 


-4.725 156 7 


-4.725 156 6 


4 


-4.630 025 


-4.630 025 


5 


-4.586 389 


-4.586 384 



TABLE IV: Comparison of the excited state energies (with errors in the least significant digit) 
of Ar4 in D = 3 and 5 dimensions. 



k 


D = 4 


D = 6 


1 


-5.821 21 


-5.821 21 


2 


-5.346 6 


-5.337 2 


3 


-5.26 


-5.18 


4 


-5.06 


-4.99 


5 


-4.95 


-4.91 



TABLE V: Comparison of ground and excited state energies (with uncontrolled errors) of Nes 
in D = 4 and 6 dimensions. 



10 



III. EXACT RESULTS 



A. Clusters in arbitrary number of dimension 



The Schrodinger equation for an N particle cluster in D spatial dimensions is a differential 
equation in ND variables. For a S state the wave function is invariant under rotations and 
translations. Therefore, one can write the wave function as a function of fewer than ND 
variables. To accomplish this we proceed as follows. 

Consider the N — 1 difference vectors r 2 i, r 3 i, . . . , r^i as defined by Eq. (2a). Note that 
these vectors cannot be linearly independent unless D > N — 1, in which case they define 
a parallelepiped P. Precisely |(iV — 1)N independent variables are required to define P up 
to a congruence transformation. One possible choice of such variables consists of: (1) the 
angles % between the the vectors r^i and r^i or their cosines 



with 1 < i < j < N; and (2) the lengths of the vectors r a with 1 < i < N. Alternatively, 
as independent variables one may choose the lengths of all distinct inter-particle distances 
r 'ij = r ji w ifb i 7^ j- These are the variables we shall use in this paper with the assumption, 
required for linear independence, that D > N — 1. 



with a Hamiltonian slightly more general than the one defined in Eq. (3) with a potential 
that is rotationally and translationally invariant, but not necessarily a sum of two-body 
contributions. Furthermore, the mass of each particle may be different. 

We restrict ourselves to S states and to cases in which D > N — 1 so that, as discussed in 
the previous subsection, the wave functions can be considered to be a function of independent 
inter-particle distances with 1 < i < j < N. 



B. Generalized Schrodinger equation 



We consider a D-dimensional Schrodinger equation of the form 




(15) 
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By straightforward application of the differential operator identity 



3& J 



one obtains 



where 



d d 2 
V- = V ai.j— + V g l;jk (17) 
9r u T?- dr i:j dr ik 



a=l UX ™ ,% 3 



and 



g ^ ^ ^"ij ^ij ' ^ik (19) 

With the inter-particle distances as independent variables, the Schrodinger equation as- 
sumes a form that involves: (1) a linear differential operator that explicitly depends on the 
spatial dimensionality D; and (2) a second-order differential and a potential energy operator 
that are independent of D, as is clear from Eqs. (17) and (19). 

Next, we transform the Schrodinger equation into an equation in which the second-order 
operator is unchanged, the linear operator is absent, and in which the potential is modified 
by an additional term. 15 This is accomplished as follows: 

H*0 = Etp — > H> = E(f) (20) 

with 

V> = X0, (21) 

and 

H' = -H X . (22) 

X 

The action of the operator on the right-hand side of an arbitrary function is to be evaluated 
from right to left, so that multiplying by x takes precedence over operating by H. 
This yields a special case of Eq. (3.8) of Ref. 16 



H' = V - 7^r(Si + Ti + Ui) (23) 

with 



2rrii 
i=i 1 



j,k^i lJ 
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^E^^w'l ^ (25) 



We define square matrices of order N — 1, 

= (fij9i;jkfik)j,kjti (27) 

for i = 1, . . . , iV. The matrix is the Grammian associated with the N — 1 vectors with 
j — 1, . . . , i — 1, i + 1, . . . , N. Its determinant 

to = det (Gi) (28) 

is the square of the volume of the parallelepiped defined by the vectors pointing from particle 
% to all other particles. This volume is equal to N\ times the volume of the (N — l)-simplex 
of which the N particles are the vertices, which explains why uj does not depend on i, as 
our notation indicates. 

In the Appendix we show that T { vanishes for the choice 

X = u (l - D)/ \ (29) 

while 

U^hN-lf-iN-Df]^^ (30a) 
(N - l) 2 — (N — Df ^ duo dco 



M S^'^' (3 ° b) 

Clearly, Ui depends on the spatial dimensionality via its D dependent amplitude, which 
is symmetric in D about D = N . Recalling that this derivation is valid only for values 
for D > N — 1, we find that for S states the Schrodinger equation has the same energy 
eigenvalues in D = N — 1 and D = N + 1 — and for those values of D only — while the 
eigenstates are related via 

t/j(D = N- 1) = yfoij)(D = N+ 1). (31) 

Furthermore, using the fact that gijk, defined in Eq. (19), is an inner product, one can 
rewrite the sum in the Eq. (30b) as a sum of squares. This implies that to linear order in 
perturbation theory the energy eigenvalues depend quadratically on D with a minimum at 
D = N, in agreement with our numerical estimates presented in Section II B. 
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IV. DISCUSSION 

We transformed the Schrodinger equation for the rotationally and translationally invari- 
ant states of an iV-particle cluster in D > N—l spatial dimensions into a differential equation 
in |(iV — 1)N independent variables in which the dependence on D is fully contained in an 
effective potential energy. Here D can be interpreted as a continuously varying parameter, 
as is commonly done in dimensional scaling studies. 17 In agreement with work by Gu et 
al., 2 we observed that there exists an inter-dimensional degeneracy, of an iV-body cluster 
in D = N — 1 and D = N + 1 dimensions. This degeneracy exists for all S states, i.e., for 
both the ground and excited states. Furthermore, the minimum energy was observed to be 
at the dimension D = N. 

We stumbled upon this inter-dimensional degeneracy numerically by Monte Carlo meth- 
ods, for which the generalization to arbitrary discrete dimensions is simple. In retrospect, 
knowing that this degeneracy is an exact property of the Schrodinger equation is useful be- 
cause it provides a powerful check of the validity of our computer code and of our estimates 
of systematic and statistical errors in our numerical results. 

Our results have another interesting implication. Our transformation of the Schrodinger 
equation to a differential equation in |(iV — l)iV independent variables is valid only for 
D > N—l. However, D in the resulting equation can be interpreted as a continuous 
variable, and the equation has an analytical continuation for D < N — 1 and is symmetric 
about D = N. As a consequence, the energy spectrum for D < N — 1 of the transformed 
equation is analytic in D and symmetric about D = N. Whatever might be the meaning of 
this spectrum obtained by analytic continuation of the spectrum for large values of the D 
spectrum, it cannot have anything to do with the physical spectrum of N particle clusters 
for D < N — 1. This in turn implies that there is a fundamental problem with the work by 
Gonzalez and Leal, 18 who have used 1/D expansion to estimate energy levels of Lennard- 
Jones clusters in D = 3 with iV = 3,4, . . . , 7 and 13. If such calculations could be done 
to infinite order and re-summed to yield a convergent expression valid for all D, the result 
would agree with this the analytic continuation discussed above, but not with the physics 
of clusters with more than four particles in three dimensions. 
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APPENDIX 

Without loss of generality we can restrict our discussion to the contribution to the trans- 
formed Hamiltonian H' of the kinetic energy of particle i = N. Correspondingly, we shall 
simplify our notation as follows: 

Oj = a N]i (A.l) 

9ij = 9N;ij (A. 2) 

Ti = r Ni (A. 3) 

Note in particular that n is not the distance of particle % to the origin, as suggested by 
convention and Eq. (1), but rather the distance of particle N to particle i. 
Define 

G = (W,)5=i- (A.4) 
Then the square of volume of the parallelepiped is given by the Grammian 

u = det(G). (A.5) 

Consider a symmetric s x s matrix M of the form 

(Ui(xi) if % — j 

(A. 6) 
VijyXi) Xjj Vji(^Xj] Xij if % <C J 

Since only row and column i depend on Xi this implies that 

= ^ Uji ^ 5ii + 2 ftj^) (1 _ ^ (A7) 

where M^- is the cofactor of M. 
Now 

rtC 1 

w -rj>, (A.8) 

% = 1, . . . , JV — 1, where C p has elements equal to unity on the off-diagonal of row and column 
p and equal for diagonal element p, i.e., 

i p) = + ^ (a.9) 



^ = 2r, (A.IO) 
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From this we find 

dull „ , . s 

fti -- = 2. (A.11) 

This can be verified by solving this last equation for J^-^- by means of Cramer's rule, which 
indeed yields back Eq. (A. 10). Thus we find 

Edui uo , , . 

*.^ = 2- (A.12) 

i J 

from which it follows that 

X = u^/ A (A.13) 

so that 

i=l J 

which shows that the linear differential operators vanish in T N . 

To make the dependence on the spatial dimension D explicit, we write the first term on 
the right-hand side of Eq. (26) as 

j=i j j=i j j 

To calculate the second term we use the law of cosines in the form 

ij * ik jk / * -, n \ 

9i;jk = -, (A. 16) 



which yields 

dgi-jk _ J. 1_ 

By repeated use of Eq. (A. 14) we find 



9i;jk (A.l 7) 



J,k=l J J,k=l J J,k=l x 



JV-1 
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